function [ W1 H1 ] = GetStartingWHProport( duomenys, R )
    
    data = duomenys;
    
    rows = size(data, 1);
    cols = size(data, 2);
    
    W1 = mean(duomenys, 2) ;
    W1 = repmat(W1, 1, R);
    
    H1 = mean(duomenys, 1) ;
    H1 = repmat(H1, R, 1);
    
    
    
    
    scale = 1;
    
    Wr = rand(rows, R) * scale;
    Hr = rand(R, cols)  * scale;
    
    W1 = W1 + Wr;
    H1 = H1 + Hr;
    
    
    
    
    ee = 0.01;
    
    W1(W1 <= 0) = ee;
    H1(H1 <= 0) = ee;
    
    
    
    W1 = W1 / max(max(W1));
    H1 = H1 / max(max(H1));
    
    
end

